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Abstract 



We construct the quantum versions of the monodromy matrices of KdV theory. The 
traces of these quantum monodromy matrices, which will be called as "T-operators" , act 
in highest weight Virasoro modules. The T-operators depend on the spectral parameter A 
and their expansion around A = oo generates an infinite set of commuting Hamiltonians of 
the quantum KdV system. The T-operators can be viewed as the continuous field theory 
versions of the commuting transfer-matrices of integrable lattice theory. In particular, 
we show that for the values c = 1 — 3 ^n+l > n = 1)2,3... of the Virasoro central 
charge the eigenvalues of the T-operators satisfy a closed system of functional equations 
sufficient for determining the spectrum. For the ground-state eigenvalue these functional 
equations are equivalent to those of massless Thermodynamic Bethe Ansatz for the minimal 
conformal field theory M.2,2n+3', m general they provide a way to generalize the technique 
of Thermodynamic Bethe Ansatz to the excited states. We discuss a generalization of our 
approach to the cases of massive field theories obtained by perturbing these Conformal 
Field Theories with the operator $1,3. The relation of these T-operators to the boundary 
states is also briefly described. 
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The studies of the last decade revealed a deep relation between the structures of 
Conformal Field Theory (CFT) 0, Integrable Field Theory @, § and Solvable Lattice 
Models H. The conformal symmetry of CFT is generated by its energy-momentum tensor 
T(u), whose mode expansion 

+00 

r («) = "^ + !>-» C<BU (!) 
—00 

is expressed in terms of the operators L n satisfying the commutation relations 

[L n , L m ] = (n-m) L n+m + j^(n 3 - n) 5 n+mj0 (2) 

of Virasoro algebra Vir p|. Here c is the central charge which is the most important 
characteristic of CFT; in writing (2) we have chosen the periodic boundary conditions 
T(u + 2%) = T(u). The "crural" space of states H c hirai is built up from a collection of 
suitably chosen irreducible highest weight modules , 

1~(-chiral = ©a V a > (3) 

where V a = Va q , the parameters A are the highest weights, and associated highest weight 
vectors | A) G V a satisfy the equations 

L n | A) = for n > ; L \ A) = A | A) . (4) 

Along with this "conventional" characterization of CFT a quite different description 
of CFT in terms of "massless S- matrix" has been recently proposed [Q, 0, 0. In this 
approach the chiral states of CFT are described as the scattering states of a collection of 
"massless particles" and the CFT is characterized 2 by their factorizable S-matrix, rather 
then in terms of Virasoro algebra and its highest weight representations. So far the relation 
between these two description is understood to a very limited extent. The "massless S- 
matrix" allows one to compute, through the Thermodynamic Bethe Ansatz technique 0, 

1 Of course the physical space of states 7i p hy S is embedded (by diagonal or some other suitable 
embedding) into a tensor product Tichirai <8> Hchirai, see e.g. ]IJ. 

2 Unlike the case of massive field theory the "factorizable S-matrix" description of CFT in 
general is not unique; there can exist several different choices || of the "massless particle" basis 
in Tichirai, with different particle contents and scattering amplitudes, which correspond to massless 
limits of different integrable perturbations of this CFT . 
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T0[1 , the asymptotic density of states and thus gives the central charge c (and sometimes 
few of the highest weights A ||11|| ). In general case even the correspondence between the 
"massless particle states" and vectors in (|3|) is not known. At the same time gaining more 
understanding about this relation is very important as the "massless S-matrix" description 
seems to capture much of the "integrable structure" of CFT, i.e. it is close to its description 
in terms of action-angle variables. 

From the algebraic point of view these "massless particle states" are nothing else 
then the eigenstates of an infinite set of commuting Integrals of Motion (IM) . It has been 
known for some while that if we consider the algebra UVir generated by the energy- 
momentum tensor ([!]) along with various "composite fields" built as a powers of T(u) and 
its derivatives, this algebra contains an infinite-dimensional abilean subalgebra fl2Tl , [jT3| 
spanned by the "local IM" /2/c-i £ UVir, k = 1, 2, ... which have the form 

hk-i = [ 7T T 2k (u) , (5) 
Jo 27r 

where the densities T2k(u) are appropriately regularized polynomials in T{u) and its deriva- 
tives. The first few densities T2k{u) can be written as 

T 2 (u) = T{u) , T 4 (u) =: T 2 {u) : , T 6 (u) =: T 3 (u) : +^ : (7») 2 : , ... (6) 

Here the prime stands for the derivative and : : denotes appropriately regularized products 
of the fields, for example 

where the symbol T denotes the "chronological ordering" , i.e. 

T(A(w)B(u)) = \i^ B J U \> !! f^u>^nw; 
K y ' K n \ B(u)A(w) , if Sm w > 3m u . v ' 

In writing (^) we disregarded all the terms which are total derivatives and do not contribute 
to (|5p. Although a general expression for all densities T2k{u) is not known they are 
uniquely determined (up to a normalization which we will fix later) by the requirement of 
the commutativity 

[l2fc-i,Jaj-i] = (9) 
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and the "spin assignment 
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<[pL(w- u)T(T(w)T 2k (u)) = 2k T 2k (u) . (10) 
Jc 2m 

More representatives (beyond (||) ) of this infinite set of densities T 2 fc(w) can be found in 
T2| . The integrals (|J) define operators l2k-i '■ Va — > Va which can be expressed in terms 



of the Virasoro generators L n , for example 



T -2 Vf T + 2 T i c ( 5c + 22 ) 

Is -2 2^L_ n L n + L - — L + 28g0 , 

71=1 

/ 5 = Yl ■ L ni L n2 L n3 ■ + ra 2 -l-|) L_ n L n + (11) 

"■i+'^2+'"3=0 n=l 

3^ c + 4 2 (c+2)(3c + 20) c (3c + 14) (7c + 68) 

r=l 



The "normal ordering" : : in these formulas means that the operators L n with the bigger 
n are placed to the right. Note that although these operators are not polynomial in L n 
their actions in Va are well defined. 

In this work we study the problem of simultaneous diagonalization of the operators 
Izk-i m Va by using the approach which can be regarded as a version of the Quantum 
Inverse Scattering Method 0, ||. It was remarked many times |0|, [JIB]] , [ |14"f that this 
problem can be thought of as the quantum version of the KdV problem as it reduces 
to the classical KdV problem (with periodic boundary conditions) in its "classical limit" 
c — > — oo. Indeed, in this limit the substitution 

Tin) - - | U(u) , [ , ] - ^ { , } 
{U{u + 2%) = U{u)) reduce the algebra ([!]), (g) to the Poisson bracket algebra 



3 Another way to formulate this requirement is to assign the grade 2 to the field T(u) and the 
grade 1 to the derivative; then the density T2k(u) is a homogeneous polynomial in T(u) and its 
derivatives of the total grade 2k. 
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{U(u), U(v)} = 2 (U(u) + U(v)) 8'{u -v) + 5"'{u - v) , (12) 
which is well known to describe the second Hamiltonian structure of the KdV equation 

(cl) 

provided we take one of the infinite set of classical IM I^k-i 



as the Hamiltonian. The classical IM (|T^) , which form a commutative Poisson bracket 
algebra {^fe-i) ^2i-i) = ^' evidently are the classical versions of the operators (|J), ([ll])- 
It is also well known |T5[] that the KdV equations 

d t2h _ l U(t 1 ,t 3 ,...) = {l2k-lMtl,t3,..)} , h=U (14) 

describe isospectral deformations of the second order differential operator 

L = d 2 u + U{u) - A 2 . (15) 

In particular, if we define the 2x2 Monodromy Matrix M(A), which belongs to the group 
SL{2), as 

(^i(« + 27r),V 2 (« + 27r)) = (ip^u) , ^ 2 {u)) M(A) , (16) 

where ipi(u), ^a(w) are two linearly independent solutions to the equation Lip = 0, then the 
eigenvalues of M(A) are involutive (with respect to the Poisson structure (|T2])) Integrals 
of Motion of the KdV flows, and the trace 

T(A) = tr M(A) (17) 

can be thought of as the generating function for the local IM (^3|) as it expands in the 
asymptotic series 



2tt 



here c x = \ , c n = (2 L„ 3 ? !! , n > 1 . 



h [T(A)] - A - J2 c n & A 1_2r 



(18) 



71=1 



2 ' ^ n 2™n! ' 

It is more convenient for our purposes to start with the first order differential operator 



£ = d u - 0'(w) H-X(E + F) , 
where 0(it) are the canonical variables with the Poisson brackets 

{0(«),0(v)} = - e(u-v) ; 
e(w) = n for 27m < it < 27r(n + 1) ; n G Z 
which are related to U (u) by Miura transform [|16| 



and E, F, H are the generators of the Lie algebra si (2), 



(19) 



(20) 



(21) 



[H, E] = 2E, [H, F] = -2F, [E, F] = 2H . 
In general, the classical field <p(u) has to be taken quasiperiodic 



(22) 



[u + 2tt) = 4>{u) + 2mp 



(23) 



to guarantee the periodicity of U(u). In order to define the monodromy matrix for the 
operator ( |i~9|) one has to pick some matrix representation for the si (2) algebra (|22|). Let 
icj [E], TTj[F], 7Tj[H] be (2j + 1) x (2j + 1) matrices, j = 0, 1/2, 1, 3/2, representing (H), 
such that 7Tj[H] = diag(2j, 2j — 1, —2j). Then, for given j, the solution to the equation 
C9(u) = is 



= 7Tj p (u)H Pexp(A / dv{e~ 2 ^ v) E + e 2 ^ v) F)) 



(24) 



where the symbol P denotes the "path ordered" exponential and is arbitrary vector in 
C • 7+ . The associated monodromy matrices have the form 



2tt 



M j (A) = tt j e 2 ™ pH Pexp(A / dv(e~ 2 ^ v) E + e 2 ^ v) F)) 



'o 
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(25) 



Let us introduce auxiliary matrices 



(26) 



They satisfy the "r- matrix" Poisson bracket algebra [17 



{Lj(A) ®Lj>(fi)} = [r jr (\fi 1 ),Lj(A) ® Tjj'(fi)] , 
where rjj>(\) = nj <g> 7Tf[ r ] is the "classical r-matrix" 

It follows immediately that the quantities 



(27) 



(28) 



Tj(X) =tr Mj(X) 
are in involution with respect to the Poisson bracket (|T2|), (p0|): 



(29) 



{T i (A),T i ,(/*)} = 0. 



(30) 



In particular, Ti(A) coincides with (|T^). 

After this brief review of known classical results we turn to the quantum case. Quan- 



tum version of the Miura transform (21) is the Feigin-Fuchs "free field representation" of 



Vir 0, [[Lj 



-P 2 T{u) =: ¥ /(«) 2 :+(l-/3V'(«) + 



24 



in terms of the free field operator 



(p(u) = iQ + iPu + — -e 



where the mode operators Q, P and a n satisfy the Heisenberg algebra 



(31) 



(32) 



[Q,P] = -P 2 ; 



[d-n, dm] — — P 2 ^n+m,0 



(33) 



(P and Q commute with a n ) and the parameter j3 is related to the central charge c as 

7 



1-c 25 -c . t . 

" = ^-V^i-- (34) 



Let T v ("Fock space") be the highest weight module over the Heisenberg algebra (|3 
with the highest weight vector | p) ( "vacuum" ) defined by the equations 

P | p) = p | p); a n | p) = for n > . (35) 

The Eq. (|31~1 ) defines the action of Vir in T v . For generic c and p the space T v is isomorphic 
to the highest weight Virasoro module Va with ]18| 



The space JF p naturally splits into the sum of finite-dimensional "level subspaces" 

F P = ®!^o4 l) ! L o ^ = (A + Jf . (37) 
The "normal ordering" suitable for this representation is implied in (^T[), : : means that 
the operators a n with the bigger n are placed to the right. Of course the IM (Rl), (TTJ) 
can be expressed in terms of the Feigin-Fuchs free field (f'(u). We adopt the following 
normalization of the operators l2k-i- 



(-l) k (3- 2k / 2 ^( : M-)r :+••), 08) 

where the omitted terms contain higher derivatives of (f(u). These operators act invariantly 
in the level subspaces . Therefore diagonalization of Iik-\ in a given level subspace 
reduces to a finite algebraic problem which however rapidly becomes very complex for 
higher levels. Here we list only the vacuum eigenvalues for the first few I^k-i, 



'2fc-l 



J-2k-l 

where 

J? oc (A) 

I™ C (A) = A 2 - A -f 



p) = 4l a _ c i(A) \ p) , (39) 



A- 


c 




" 24 ' 


A 2 


c + 2 




12 


A 3 


c + 4 
8 


A 4 


c + 6 




6 



2880 

A2 (c + 2) (3c + 20) _ c(3c + 14) (7c + 68) 

576 290304 ' [ } 

/r c (A) = A 4 - — A 3 + 15c 2 + 194c + 568 
6 1440 
(c + 2) (c + 10) (3c + 28) A c (3c + 46) (25c 2 + 426c + 1400) 



10368 24883200 
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After these preparations we can define a quantum version of the Monodromy Matrices 
( p5| ) and the operators (p6|), (p9|) ; we are going to use the same symbols Lj, Tj for the 
quantum counterparts of (^6|), (|29|). Consider the following operator valued matrices [^0 



Lj(A) = 7Tj 



e™ PH Vexp(\ / e~ 2 ^ u) : : e W : q~%F)) 

Jo 



(41) 



where the vertex operators 



: e ±2 ^ u) : 



p±f3 2 



(42) 



are defined as 



: e ±M«) : = exp ( ±2 V" ^=V™) e xp( ± 2z (Q + Pu)) exp( T 2 V ^e" 1 



(43) 



n=l 



n=l 



and -E, F and iif are the generating elements of the quantum universal enveloping algebra 
U q (al(2)) 



[H,E]=2E, [H,F] = -2F, [E,F] 



q H _ g -H 
q-q- 1 



(44) 



with 



q = e 



i-n-0 2 



(45) 



The symbol 7Tj in (^TJ) stands again for the (2j + l) dimensional representation of U q (sl(2)) 
so that (HTI) is in fact (2j + 1) x (2j + 1) matrix whose elements are the operators in 



— ®n=-oo^ p+nf3 2 ■ 

These operators are understood as the power series in A, 



(46) 



Iy(A) =7T j 



iirPH 



k=0 



K(ui)K(u2)...K(uk) duidu2---duk 



1-K~>U\>U^ > . . . >ltl, >0 



(47) 



K(u) =: e~ 2 ^ u) : q% E+ : e 2ip{u) : q~^F 
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The integrals in (H7]) are convergent for 



-oo < c < -2 . (48) 

In fact, the operators (|4T| ) can be defined for wider range of c by appropriate regularization 
of divergent integrals in (|4~^) . In this paper we restrict our attention to the domain fl4~8|) . 

The operator matrices (|4j) are designed in such a way that they satisfy the Quantum 
Yang-Baxter Equation 



R jy (A^ _1 ) (Lj-(A) <8> l) (1®Ljv(/x)) = (l®L f ([j,)) (L,(A) <g> l) R^A// -1 ; 



(49) 



where the matrix R.,y(A) is trigonometric solution of the Yang-Baxter equation which 
acts in the space iXj ® ixy . In particular 



Rii(A) 



(q~ x X - qX~ x 



\ 



A -A" 1 q~ x — q 
q~ l -q A -A" 1 



\ 



q^X-qX' 1 / 



(50) 



The validity of (|4^) can be checked explicitly for the first few orders of the expansion of 
( fl9D in A and fi with the use of (flTD. One can prove (^9|) to all orders by taking the 
discrete approximations to the P-ordered integral in ([4~1~1) pOU . 
Let us define now the operators 



T J (X)=tr 7rj (e^ PH L,(A)) 



(51) 



which satisfy 



[T j (X),T j ,(n)} = 



(52) 



as a simple consequence of (^9|) . It is easy to see that the operators ( pT|) commute with 
the operator P and hence they act invariantly in T v . Moreover, one can check by direct 
calculations that 



[T i (A),J 2fc _i]=0. 

It follows, in particular, that the level subspaces Tp are the eigenspaces of Tj(A). 
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(53) 



Let us concentrate first on the simplest nontrivial T-operator T(A) = Ti(A) which 
corresponds to the two-dimensional representation of U q (sl(2)). In this case 



Substituting (f47|) into ([51]) and computing the trace one finds 





1 



(54) 



T(A) = 2 cos(2ttP) + x2n Qr 



n=l 



where 



(55) 



Q n = q n I (e 

' 2lT ~>Ul ~>U 2 ~>...~>U 2n >0 



2iirP . e ~2tp(u 1 ) .. e 2ip(u 2 ) .. e -2ip(u 3 ) . . ^2ip{u 2n ) . _|_ 



e -2i7rP . e 2(p(ui) .. e -2ifi(u 2 ) .. e 2ip(u z ) . ^ . g - 



■^(u 2n ) . j dui ^ du 



2n ■ 



(56) 

So, the operator T(A) can be considered as the generating function for the "nonlocal IM" 
Q n which commute among themselves 



and also commute with all the "local IM" ^k-i 



(57) 



[l2k-l,Qn] =0 . 



(58) 



The operators Q n invariantly act on each of the level subspaces Tp ; in particular, the 
vacuum state | p) is the eigenstate of all Q n , 



Qn\p) = Q { r c \p)\p) 



(59) 



where the eigenvalues Qn ac \p) are given by the integrals 



i>2iT /-Mi rVi r-u 2 l>v n -i ru n 

Qn ac) (p) = dm dvi / du 2 / dv 2 ... du n I dv r 

Jo Jo Jo Jo Jo Jo 



n 

n 



2 sin( 



Ui — Uj 



Y[h s in(^— ^)) 2 cos(2p (n + ^ivi-Ui))) . 



2/T 



-2/3 



o • f v i- v 3\ 
2 sm (— r-^-) 



2/3 J 



n 

n(2-»(^) 



-2/3^ 



(60) 
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In particular 



(v) ~ 4 " 2 ^ ~ 2 ^ (61) 



Using the power series expansion (^3]) one can show that the operator T(A) is entire 
function of A 2 (just as it was in the classical case) in the sense that all its matrix elements 
and its eigenvalues t(A) are entire functions of this variable. In fact, the eigenvalues £(A) 
exhibit essential singularity at the infinity as the result of accumulation of zeroes of these 
functions along the real axis in A 2 -plane as A 2 — > — oo 4 . Asymptotic form of the operator 
T(A) at A 2 —>■ oo is of primary interest because, in view of (|18D, it is in this limit one 
anticipates to get in touch with the "local IM" (||), fl38|). Rough estimate of the matrix 
elements of the operators ( |56|) along the lines proposed in [^2| as well as physical arguments 



based on the thermodynamic treatment of the "Coulomb gas partition function" fl60| ) gives 
the leading asymptotic 

log[T(A)] ~ m A 1+c , (62) 

where 



e=T^" (63) 



A 2 

1-/3 S 

The constant m is expected to depend on c but not on the particular matrix element of 
this operator; this asymptotic form holds in the domain 



-7r + e < arg( A 2 ) < n - e, A — ► oo (64) 

with arbitrary small positive e. So far we could not find any direct way to compute 
the constant m. Therefore we used the Algebraic Bethe Ansatz to study the discrete 
approximations to the integral in (f4~T|); in taking the continuous limit the integral form of 



the Bethe Ansatz equations proposed in [23[ is particularly useful. This way we obtain 



T(A)=A( ? A)+A- 1 ( ? - 1 A) , (65) 
where A (A) expands in the domain fl64|) into the asymptotic series 

oo 

logA(gA) ^ m A 1+ « - ^ C n J 2n _j X^- 2n ^ 1+ ^ , (66) 

n=l 



Validity of this picture is restricted to the domain (| 
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where 




(67) 



C " ~" "\ \ l + U \ in T(l-l)l r(l + (n-|)e) 



and l2k-i are exactly the "local IM" (||), ([38D. The coefficients C n here are closely related 
to the canonical on-shell normalization of the local IM of the quantum Sine-Gordon theory 
|24{ j. Note that the coefficients c n in (|T3|) can be recovered in the classical limit £ — ► 
(c — > — oo). Appearance of the fractional powers of A in fl66| ) should not be very surprising 
- the exponentials in (^1]) have anomalous dimensions and so the spectral parameter A is 
to be thought of as carrying the dimension [ length ]~~- The expansions (|55]) and (|66| ) 
give highly nontrivial analytic relation between the "nonlocal IM" (|5B) and "local IM" 



The higher spin operators Tj (A) also admit the power series expansion similar to (|55|), 
the coefficients being of course algebraically dependent from the operators Q n in ( |55|) . 
For the first few coefficients we find 

, . sin (2ttP (2j + 1)) 9 A . 9 , 

sm (2ttP) (68) 

A 4 [ Aj(2irP, 2n/3 2 ) Q 2 + B j (2ttP, tt/3 2 ) Q\ ] + 0(A 6 ) , 
where numerical coefficients read explicitly 

1 / sin(2j + 1) (x - a) sin(2j + 1) (x + a) 



Aj(x, a) 



4 sin x sin a I sm(x — a) sin(x + a) 



Bj (x, a) = : I ^(2j + l)(x-2a) 

16 sin x sin a sin 2a \ sin(x — a) sin(x — 2a) 



sin(2j + l)(x + 2a) 2 sin(2j + l)x cos a 



sin(x + a) sin(x + 2a) sin(x — a) sin(x + a) 

The operators Qi and Q2 in are the same as in (|56|) . This algebraic dependence can be 
summarized by the functional relations 

T^A^g-U) = l + T j _ h (\)T j+ i(\) . (69) 
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Note that these relations are identical to the functional relation obeyed by the commuting 
transfer-matrices of the integrable XX Z model ||25|| , ||26|| , |[27|| . This is well expected as 



our T- operators appear to be the continuous field theory versions of the lattice transfer- 
matrices. The relations ( |69|) can be derived from (|49| ) essentially the same way they are 
obtained in the lattice theory, by using the R-matrix fusion procedure ||21|| . The simple 



form of flol^D is due to the fact that quantum determinant |J of the operator L 1 (A) ( ]4T] ) 
is equal to 1. Using the relation (^9[) one can show that all the operators Tj(A) with j > 1 
are also entire functions of A 2 as well as T(A). 

At generic values of c the relations (^9|) just allow one to express the higher T 
operators in terms of the lower ones. However it is well known in the lattice theory 
that at particular values of the parameters, when the XX Z system can be reduced to the 
RSOS model p8|, the functional relations truncate to become a finite system of functional 



equations |2"5 |, Of course, similar phenomenon happens in our continuous theory. 

Restricting our attention to the domain (|4^) we find that the most simple truncation 
occurs at 

(2n + l) 2 2 
c=l-S \ ' , ^ = ^T— ?! n = l,2,3,.... (70) 
2n + 3 2n + 1 

For given n in (|7(]) consider the finite collection of the Fock spaces 

2/c-2n-3 , _ 

^ ! Pk = 2(2n + 3) - k = l,...,n+l. (71) 



At these values of p the Fock spaces ([71]) are known to be reducible with respect to the 
action of the Virasoro algebra (they correspond to the (1, k) degenerate representations in 
the Kac classification [ffQ| , Ifl8| ). Let us denote V Pk the associated irreducible Virasoro 
module obtained from (|7T| ) by factoring out all the submodules. Then the space 



coincides with the space of chiral states of the "minimal CFT" A^2,2n+3- R is possible to 
show that the operators Tj(A) with j = 0, |, 1, n+ | invariantly act in the space (|7^ ) 
and being restricted to this space these operators satisfy the symmetry relation 



T n+ i_ i (A) = T i (A); i = 0,i,..,n+i , 
14 



(73) 



in particular T n+ i (A) = I. Under these circumstances the relations (^9[) become a finite 
system of functional equations 

tj(q^) tj(q-*\) = l+t j+ i(X) tji(\) ; 

2 (74) 
to (A) = (A) = 1 ; Wl-^A) = t;(A) , 

which is satisfied by all eigenvalues tj(X) of the operators Tj(X) in the space (|7^). We 
conjecture that any solution t,(A) to the functional equations ( |7^ ) which is an entire 
function of A 2 and has the asymptotic behavior (^2|) in the domain ( p^ ) corresponds to 
some eigenstate of Tj(A) in the space Q7^). 

It was recognized before [^7] that the substitution 



Y j (9)=t jH (X) ^_i(A) ; A = exp( r ^) 
brings the system ( |7^ ) to the form 



(75) 



Y, (0 + ^) Y 3 (e-^) = (l + Y jH (9)) (1 + Y s _ k {ft) ; ^ 
r o (0) = r„+i (*) = o ; = , 

which coincides with the functional form of the Thermodynamic Bethe Ansatz (TB A) equa- 
tions lyj for the "massless S-matrix" theory associated with the minimal CFT M.2,2k+i 



|3^| . In this "massless TBA" approach one assumes that the states of the CFT M.2,2k+i 
in infinite volume can be interpreted as the scattering states of a collection of n mass- 
less right-moving particles Aj j = \ , 1 , | , . . . , ^ , with the energy-momentum spectrum 
E = P. Parameterizing the energy-momentum for the particle Aj in terms of rapidity 6 

as 

EM = PM = ^e\ mj = ^ cot(^) sin « j) • (77) 

one conjectures the purely elastic and factorizable S-matrix for these particles with the 
two-particle elements Sjj'(9 — 9') (describing the scattering processes Aj(6) + Aj/(9') — > 
Aj{9) + Aj>(9'), in obvious notations) given by fg3|, fS|] 

2min(j,j')—l 

Sjf{9) = F J+r (9)F lj _ j/l (9) n F^_ rl+k (9) , (78) 

k=i 
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where the notation 

p sinh6> + i sin(7r£ x) 

sinh# — i sin(7r£ x) 

is used. The S-matrix (|TS| ) allows one to determine the spectral density of states in 
this infinite-volume system and hence to compute the free energy of this system at finite 
temperature R . As the thermal ensemble corresponds to the circular compactification of 
"imaginary time" one can interchange the roles of the space and the "imaginary time" and 
reinterpret this free energy as the ground-state energy eo(R) of the finite- volume system 
with periodic boundary condition defined on the circle with the circumference R 5 . This 
way one obtains 

= -E ^ fj~^ d log(l + e-«<'>), (79) 
where the functions 6j(9) solve the integral equations of TBA: 

Rm i e o = ej{ e) + J2 f ^nAe-V) log(l + e-^')) . (80) 



2 

r 



with the kernel 



(fjji(6) = —i dg In Sjj/ (9) 



By inverting the integral operator in the right hand side of (|80f) one can show, following 
113 111, that the functions 



Yf- St \6)=e e ^ e) (81) 

satisfy the functional equations ( ff6|) and agree with the asymptotic condition ( |62D pro- 
vided the parameter m in (|77|) is the same as in (^2]). The ground state (the state of the 
lowest eigenvalue of L ) of the CFT A^2,2n+3 is the vacuum state | p n +i) from ([72]). We 
have checked (for the case of ^2,5) that for this value of p the numerical solution of the 
integral equations (|80|) (which is named the "kink solution" in jl0|) matches perfectly 
both the expansions (^6|) and (|68| ) [p5| . We see that the functions €j(9) of TBA admit 
an interpretation in terms of the eigenvalues of the T -operators in the "cross-channel" . 



5 Of course R can be arbitrarily changed by the scale transformation which is the symmetry 
of CFT. Note that in our previous discussion of the T-operators this parameter was set equal to 
2tt. 
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Let us stress that while the solution to the integral equation ( |gOD is unique, the func 



tional equations ( |7q) admit infinitely many solutions even in the class of entire functions 
with the asymptotic behavior (|62| ) . From the analytic point of view these solutions (which 
are in one-to-one correspondence with the eigenstates of the T operators in the space ( |72|) , 
by our conjecture) can be characterized by the patterns of zeroes of the functions tj(X) 
in the complex A 2 plane. All the solutions have infinitely many zeroes at the negative 
segment of the real axis in this plane (with — oo as their accumulation point) and some 
finite number of zeroes away from this locus. It is characteristic feature of the ground state 
solution ( |8lD that all zeroes of the functions t 9 - r ' st {\) are located at SmA 2 = 0; 5feA 2 < 0; 
under this condition flS0| ) follow from (|76|). The eigenvalue functions tj(X) associated 
with the other eigenstates in ([72J) have more complex patterns of zeroes. For example, for 
the vacuum states | pk); k = 1, 2, ...,n the eigenvalue functions tj(X) have finite number of 
zeroes at QmX 2 = 0; SfteA 2 > and no complex zeroes, and for the excited states in (|7^ ) 
these functions exhibit also complex zeroes which can be interpreted as the rapidities of 
the massless particles in these states 6 . In all cases one can use the technique developed in 
ETJ] to convert the functional equations ([7Bp into integral equations similar to ( |gOD but 
with additional terms in the left-hand sides. We have analyzed numerically these "excited 
states TBA equations" for few simplest states in the CFT and found again a perfect 



agreement with the expansions (|66|) and (^8|) |35 



Strictly speaking, the above discussion concerned the case of CFT and associated 
massless TBA. It is known however that the integrable structure of CFT remains essentially 
intact in more general massive Quantum Field Theory obtained by perturbing this CFT 
with the relevant primary field $i 5 3 ||37|| . Namely, this perturbed field theory exhibits two 
infinite sets of commuting local Integrals of Motion, l2k-i and /2fc-i, k = 1,2, where 
/2/c-i are obtained by appropriate deformations of the "right" conformal IM (|5|) and /2/c-i 
are the deformations of the same IM from the "left" sector of the CFT (see ||37|1 ). There 
are good reasons to believe that the above operators Tj(A), with suitable deformations 7 , 



It would be extremely interesting to find out how the classification of states by the patterns 



of zeroes relates to the "fermionic sum representations" of the Virasoro characters [36]. 

7 In fact, it is not very difficult to figure out what these T-operators in the perturbed theory 
should look like, if one recalls the known relation between the integrable structures of KdV and 
sine-Gordon theories [17|, [12|, [13|, [14| and the relation between the quantum Sine-Gordon theory 
and the perturbed minimal CFT explained in [38], [ pi) . As we did not yet elaborate all the details 
we do not present these deformed T-operators here. 
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can be extended to the case of the perturbed theory in such a way that they satisfy the 
equations (|5^) and (|69|) . They also enjoy the same asymptotic expansion (|66|) in terms 
of (deformed) "right" IM I^k-i- The major difference is that the deformed operators 
Tj(A) are no longer entire functions of A; instead they are expected to have an essential 
singularity at A = (being regular everywhere else in the finite part of the A plane) which 
is controlled by the asymptotic expansion in terms of the "left" IM I^k-ii similar to (|66|), 
with A replaced by fi/X, where /i is related in a simple way to the coupling constant of the 
perturbed theory In the case of perturbed minimal CFT A^2,2n+3 the functional relations 
( |69D become again the closed set of functional equations flT4] ) , but now one has to look for 
the solutions which have essential singularity at A = and enjoy the asymptotic behavior 

rri / \ \ f rn A 1+? , if A — > oo ; . . 

^ T ^)~{ m(s) ^ ifA ^ . (82) 

Like in the case of CFT, it is very plausible that there is one to one correspondence 
between the solutions of ([74]) with these analytic characteristics and the stationary states 
of the perturbed CFT A^2,2n+3 on finite circle. We studied numerically [|35| few simplest 



solutions in the case of perturbed A^2,5, again with the excellent agreement with the data 



available through the Truncated Conformal Space method ||39|| . 

We can not resist the temptation to mention here a remarkable relation between the 
operators ( |5l|) and CFT with non-conformal interactions at the boundary. According to 



the analysis in |RJ , in a minimal CFT (the one where the sum (|3|) contains finite number 
of terms) with a boundary only finitely many conformally invariant boundary conditions 
(CBC) is possible. These CBC are in one to one correspondence with the terms in (|3|); 
we denote B a the CBC associated with V a . The boundary state | B a ) e TC p h ys describing 



CBC B a has the form 



b v Ok fj. 



where S is the matrix of modular transformation of characters, {| a M )} is an arbitrary 
basis in V a and {| a^)} is the dual basis in V a . Furthermore, with a given conformal 
boundary condition B a , the space of local boundary operators 8 forms the subspace H-B a 



We call here local the boundary operators which exist as the insertions at the points where 
two components of the boundary with the same boundary conditions meet, in contrast with the 



operators associated with the juxtapositions of different boundary conditions, see [40] . 
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in Tichirai isomorphic to the direct sum ©^V;, where only those b which have nonvanishing 
fusion constants N% a are admitted. One can obtain more general non-conformal boundary 
conditions by perturbing the CBC B a with relevant operators ipf, G Ti.B a - Let us denote 
B a (g) the CBC B a perturbed with the operator -01,3, the primary field associated with the 
highest weight vector in Vi^; here g stands for the coupling parameter. The corresponding 
boundary state is 



p2n 

B a{g)) = Pexp( - g du ^1,3 (w)) I B a 
Jo 



> 



(84) 



where we have assumed the geometry of a half-infinite cylinder with u being the coordinate 
along the boundary, ^1,3 ( M + 27r) = ^1,3 It is known ||]J that minimal CFT with 
the perturbed boundary condition B a (g) is integrable. Now we recall |4^] the standard 
correspondence between the vectors in Ti c hirai ® Wchirai and endomorphisms of 7i c hirai 
and denote Hj(g) the operator in H c hi ra i associated with the boundary state ( j54"D with 
a = (1, 2j + 1). Using (^3|) with the explicit form of the matrix S for the minimal model 
M n ,n' ( /3 2 = -^r, n' > n = 2, 3, ...) we obtain 



Bi(o)=x; ( 



k,i 



/nn 
\~8~ 



sin ( 



Pk,i) 



sin (27T p k: i) 



sin (2tt p ki i (2j + 1)) P M 



(85) 



where P& / are the projectors on the subspaces V k ,i and 



2 p k ,i =l3 2 l-k 

are the eigenvalues of the operator P on these subspaces. Comparing this expression with 
IBD we see that 



3,(0) 



nn 



sm(2nP) sin(27r/r 2 P)) 2 Tj(0) 



(86) 



Of course this relation is not very surprising as the operators Tj(0) are essentially the 
Verlinde operators (p a [42] for a = (1, 2j + 1) and the eigenvalues of the latter are known 



to be Sab/Sob- What is less trivial is that the P-exponential in Q4"7| ) gives exactly the 
Feigin-Fuchs realization for the P-exponential in (|84j) and so this relation extends to the 
case of full boundary state TSj(g), i.e. 



7171 



sin (2ttP) sin (2vr/r 2 P) ) T^A) , 



(87) 
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where the relation between A and g depends on the normalization of the boundary field 
"01,3 ( u ) m (HD- Using the operator product expansion 



tfji^(u)ifji^(u') = (u — u') 2Al > 3 + less singular terms 
to fix this normalization we find 
9 2 = X 



2 _ 4 sin (2n(j + l)/3 2 ) sin (2njP 2 ) r 3 (l - (3 2 )T(3(3 2 - 1) 



7rsin(27r/3 2 ) 1-2/3 2 

The relation fl87|) allows one to interpret many properties of the T-operators discussed 
in this paper in terms of Renormalization Group flows between different CBC. We believe 
it also throws new light on the functional equations for the boundary partition functions 
obtained in |2l| by TBA approach. More details about this relation will be given in [35]. 



The above discussion was restricted to the domain (f48[) of c as for c > —2 the integrals 
in (^7|) become divergent. However this limitation is not very significant. It is possible to 
change the integration contours in (|47|) in such a way that the operator (f|l]) will make 
sense for any c (this can be thought of as analytic continuation in c). Moreover, with this 
redefinition most of the above properties of the Tj(A) operators will still hold in larger 
domain — oo < c < 1, the most significant change being in the nature of the essential 
singularity of Tj(A) at A = oo for —2 < c < 1 ( (|66|) remains valid but in the smaller 
domain — 27r(l — (3 2 ) < arg( A 2 ) < 2it{1 — (3 2 ) while completely different asymptotic 
behavior emerges in the complimentary sector). Although we did not yet completed the 
analysis for this domain —2 < c < 1, for the cases of unitary minimal models M. n ,n+x [El 



we expect to observe the "truncation" of the functional relations similar to the one 

for A-^2,271+3 discussed above. We hope to return to this point elsewhere. 

And, of course, the above approach can be generalized in a straightforward way to 
describe the quantum theory of generalized KdV [13] associated with the higher-rank 



simply laced Lie algebras, or equivalently, to describe the integrable structure of the CFT 
with the extended W-symmetry. 
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